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We derive analytical formulas for the forward emission and side emission spectra of cavity-modified 
single-photon sources, as well as the corresponding normal-mode oscillations in the cavity quantum 
electrodynamics strong-coupling regime. We investigate the effects of pure dephasing, treated in 
the phase-diffusion model based on a Wiener-Levy process, on the emission spectra and normal- 
mode oscillations. We also extend our previous calculation of quantum efficiency to include the 
pure dephasing process. All results are obtained in the Weisskopf-Wigner approximation for an 
impulse-excited emitter. We find that the spectra are broadened, the depths of the normal-mode 
oscillations are reduced and the quantum efficiency is decreased in the presence of pure dephasing. 
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I. INTRODUCTION 

Single-photon sources (SPS) have important uses in 
quantum communication [j], [2j , quantum computing Q , 
and metrology [4| . Since the turn of this century, signifi- 
cant progress has been made in generating single photons 
from a microscopic quantum emitter such as an atom or 
ion [1, @] , an organic molecule @, H, , a semiconductor 



uantum dot (QD) or nanocrystal (Tajn|JT3JTE Q EI 



16|,ll7|, or a color center in diamond^ilTll^Most of 
these SPS are based on spontaneous emission, whose life- 
time eventually limits the emission rate and linewidth of 
the SPS, and whose isotropic nature prevents high collec- 
tion efficiency. A more promising scheme for producing 
well-controlled single photons is cavity quantum electro- 
dynamics (QED). An atom or QD inside a high- finesse 
microcavity is prepared in an excited state and is allowed 
to spontaneously emit. The emission rate and linewidth 
can be considerably altered [HI, [l3[ in the cavity-QED 
weak-coupling regime. In the bad-cavity limit, the al- 
tered emission rate is obtained from Fermi's golden rule 
using a modified density of states to account for the cav- 
ity boundary conditions. In the weak-coupling regime, 
the atomic excitation is irreversibly lost to the contin- 
uum of all available photon states, including the side 
modes (leak modes) and cavity modes. The side emis- 
sion and forward emission spectra are single-peaked, ei- 
ther enhanced or inhibited. 

If the cavity volume is sufficiently small, and the finesse 
is high enough such that the coherent interaction rate be- 
tween the quantum emitter and cavity exceeds the decay 
rates of the composite system to be in the cavity-QED 
strong-coupling regime, the emission process is reversible 
and a photon emitted into the cavity can be coherently 
reabsorbed before it is emitted out of the cavity. The ini- 
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tially excited quantum emitter undergoes single-quantum 
Rabi oscillations. The emissions, still allowed to both the 
side and forward direction of the cavity, will show double- 
peaked spectra. In this regime, the coupling of the emis- 
sion to the single cavity mode, however, is far stronger 
than its coupling to the side modes. If, in addition, there 
is almost no dephasing of the quantum emitter during 
the emission process, the emission process can be nearly 
deterministic , emitting a photon into a well-defined 
"forward" beam outside the cavity. 

Most if not all single-quantum systems, however, in- 
evitably interact with certain heat baths, leading to de- 
phasing or loss of coherence, which results from a ran- 
domization of the phases of the emitter's wave functions 
by thermal fluctuations in the environmental fields. Pop- 
ulation relaxation processes contribute to dephasing with 
a dephasing rate given by half the population decay rate. 
It is often necessary to account for other dephasing in- 
teractions, such as elastic collisions in an atomic vapor, 
or elastic phonon scattering in a solid, the so-called pure 
dephasing process. Pure dephasing causes the coherent 
overlap of the upper and lower state wave functions to 
decay in time, while not affecting the state populations. 
For example, the pure dephasing rate can be small and 
ignored for resonant excitation of a single QD at low tem- 
perature (6 K) and power density [2111 . While at elevated 
temperature, however, experiments |22Ll23| reveal a pure 
dephasing contribution that dominates excitonic dephas- 
ing. Our results are directly applicable to experimental 
data presented in Refs. [Ij, [H HH, for which no theo- 
retical predictions were previously available. 

In this paper, we derive the normal-mode oscillations 
and spectra of the photon emitted both in the side di- 
rection and in the forward beam in the cavity-QED 
strong-coupling regime. We first review the derivation 
of the normal mode oscillations and emission spectra ob- 
tained in the Weisskopf-Wigner approximation (WWA) 
[13] in Sees. [Til an d IIIII Then we focus on the influ- 
ence of the pure dephasing process, treated in the phase- 
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detunings of the emitter-cavity, emitter-reservoir, and 
cavity-reservoir, a and fit are the annihilation and cre- 
ation operators for the single cavity mode under consid- 
eration, while a z and a± are the Pauli operators for the 
emitter population inversion, raising, and lowering, re- 
spectively. Here we treat the atomic transition frequency 
ujq as constant. Later we allow it to fluctuate, to model 
pure dephasing. 

Given that there is only one excitation in the system, 
the state vector can be written as 



FIG. 1: (Color online) Schematic description of a lossy two- 
level emitter interacting with a single mode in a leaky optical 
cavity, go is the coupling constant between the emitter and 
the cavity field. Ap, Ap and £>g, _B| are the coupling con- 
stants between the emitter, a single photon and their respec- 
tive reservoir fields (Ri, R2)- 
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diffusion model based on a Wiener-Levy process [28[ , on 
the normal-mode oscillations and emission spectra in Sec. 
TVl We also extend our recently presented results for 
single-photon quantum efficiency (QE) [29[ to include the 
effects of pure dephasing of the atomic dipole. We find 
that the depths of the normal-mode oscillations are re- 
duced, emission spectra are broadened and the QE is 
decreased in the presence of pure dephasing. 



II. PROBABILITY-AMPLITUDE METHOD IN 
THE WEISSKOPF-WIGNER APPROXIMATION 

Consider the interaction of a quantized radiation field 
with a two-level emitter, an atom or QD, located at an 
antinode of the field in an optical microcavity with a 
length L, as in Fig. [TJ Mi is a perfect 100%-reflecting 
mirror and M2 is a partially transparent one, from which 
a sequence of single photons-on-demand emerges. 

The interaction Hamiltonian Hi{t) in the interaction 
picture for this system in the dipole approximation and 
rotating-wave approximation is [301 ] 



Hi(t) =hg (&+ae lAt + H.c.) 



(1) 



where |m, n) (to = e, <?; n — 0, 1) denotes the emitter 
state (excited state, ground state) with n photons in the 
cavity. \jp) r^I^) r 2 (ji ^ = 0, 1) corresponds to j photons 
in the p mode (other than the privileged cavity mode) 
of the emitter reservoir R\ and I photons in a single- 
mode (k) traveling wave of the one-dimensional photon 
reservoir R 2 (output beam). E(t), C(t), Sff(t), and 0%(t) 
are the slowly varying probability amplitudes. 

The equations of motion for the probability amplitudes 
are obtained by substituting \ijj(t)) and Hi(t) into the 
Schrodinger equation and then projecting the resulting 
equations onto different states respectively. In the WWA 
1271 [3fj| , we obtain 



E(t) = -ig e* At C(t)- 1 E(t), 
C(t) = ~ig e~ lAt E(t) - nC(t) 

S fi (t) = -iA} ( dt" e lS ^" E{t") 
Jo 

O g (t) - -iB% j dt'e lSkt 'C(t') 
Jo 



(3) 



(4) 



where A = ujq — lu c , S p = lu„ — luq, 5^ = lj^ — lo c are the 



where 7 and k are one-half the radiative decay rates of 
the atomic population (other than the privileged cav- 
ity mode) and the intracavity field, respectively. Dots 
indicate time derivatives. The general solutions to the 
coupled differential Eqs. ^ and ([J} are 
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k + 7, r 



7, 



and A 



Vflg-[(r-iA)/2]2._ 

In the strong-coupling regime, defined by go 3> k j 7: 
the real part of A is much larger than its imaginary 
p art. Then A can b e approximated as A rs g = 
yjgl + (A/2) 2 - (r/2) 2 , which is the generalized vacuum 



Rabi frequency. Note that for the case when the emitter 
and cavity are exactly at resonance, A = loq — oj c = 
0, the comple x frequency A is purely real and equals 
\J g^ — (r/2) 2 . The solutions to the probability ampli- 
tudes arc then 
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In this context, we assume the quantum emitter is pre- 
pared in an excited state E(0) = 1, C(0) = at time 
to = (more generally, it can be prepared in an arbi- 
trary single-quantum state). The solutions subject to 
this initial condition are 
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Sp(t) and Oz(t) can be obtained by carrying out the in- 
tegrations in Eq. (01 . 



III. NORMAL-MODE OSCILLATIONS AND 

EMISSION SPECTRA OF SPS IN THE 
CAVITY-QED STRONG-COUPLING REGIME 

The strong interaction between an excited quantum 
emitter and a single cavity mode leads to single-quantum 
Rabi oscillation (normal-mode oscillation) in the time do- 



main or a frequency splitting in the frequency domain, 
the so-called normal-mode splitting, which arises from 
the coherent interaction of two degenerate systems — 
the single quantum emitter and the single cavity mode. 
In this section, we discuss the normal-mode oscillations 
and emission spectra of SPS in the cavity-QED strong- 
coupling regime. We first investigate the normal-mode 
oscillations by calculating the probabilities of finding the 
composite system in different states. Then we define and 
calculate the spectra of SPS appropriate for this case in 
the long-time limit. 



A. Normal- mode oscillations 

The normal-mode oscillations can be viewed either in 
the dressed-state picture or in the bare-state picture. 
Here we look at them in the bare-state picture where the 
oscillations can be relatively easier to find. The proba- 
bility of finding the system in the excited atomic state 
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The probability of finding the system in the single cavity 
mode is 

P c (t) = \C(t)\ 2 = 4e- Kt sin 2 ( 3 i). (10) 
9 

Consider the case when the emitter and cavity are exactly 
at resonance, A = ojq — ui c — 0. Figure [2] are plots of the 
two probabilities with both linear and logarithmic scales. 
The probabilities oscillate sinusoidally with an exponen- 
tial decay envelope. However, they have opposite phases, 
which indicate the coherent oscillatory energy exchange 
between the excited emitter and the cavity field. 



Define the emission probability P (t) to be the prob- 
ability of finding a single photon in the output mode of 
the cavity between the initial time to = and a later 
time t. This equals 



P {t) = 2k / dt'\C(t')\=rj q 
Jo 
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where r) q = [go/(g 2 + ivy)] [k/(k + 7)] is the single- 
photon QE, given by the single-photon emission prob- 
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FIG. 2: (Color online) Probabilities of finding the system in an excited atomic state and in a single cavity mode, red square 
and blue dot curves respectively, with (a) linear scale and (b) logarithmic scale, given (go, k, j)/2-k — (8.0, 1.6, 0.32) GHz. 



ability P (t) in the sufficiently long-time limit t ^> K 1 , 
which has been discussed in our previous publication [2^ | . 

B. Emission spectra 

It may seem strange to talk about the spectrum of 
a single-mode field since we normally associate a single 
mode with a single frequency. Here we are dealing, how- 
ever, with what should more correctly be called a quasi- 
mode, a mode defined in a leaky optical cavity, which 
therefore has a finite linewidth. For a stationary and er- 
godic process, the Wiener-Khintchine theorem [31J states 



that the spectrum is given by the Fourier transform of the 
two-time correlation function of the radiated field. In the 
strong-coupling regime and for an impulsive excitation of 
the system, however, this relation between the correlation 
function and spectrum fails because the coherent inter- 
action overwhelms the relaxations here. There is no time 
t after which the correlation functions depend only on 
the time difference. Thus the dipole correlation and the 
emitted field correlation cannot be stationary. We use a 
generalized definition of the Wiener-Khintchine spectrum 
appropriate in this case (Appendix 

The side emission and forward emission spectra are 
defined as, in the long-time limit (t 3> K~ r ) 



SseW) = V R °{/ dT£m ' T J dtE ( t + T ) E *( t ) 

dtC(t + T)C*(t) 



S FE {ty = —Fie 



dre 1 



(12) 
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where f2' = u> — luq is the side emission frequency centered 
at the atomic transition frequency luq, and £1 = lu — lo c 
is the forward emission frequency centered at the cavity 
resonance uj c . Note that W = fl — A, they are equal 
when the atom-cavity is at resonance. Substituting the 
Eqs. © and © into Eqs. flTJJ) and (US]), we obtain the 
unnormalized spectra 
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When writing them explicitly in terms of frequency 
variables u, u>o, and u> c instead of fi', O, and A, Eqs. 
dHJ) and CG>]) change to 



S S e{u - cj ) = - 



Sfe(u — w c ) = — 

7T 



k — i(u> — LU C ) 



[K/2 - i(2u) -u) - u c )/2Y + g 2 



[K/2 - i(2w - w - w c ) /2f + g 2 



The side emission and forward emission spectra are typ- 
ically measured independently in experiment. The de- 
tuning A in two cases are realized differently. For ex- 
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ample, to measure the side emission spectrum, we keep 
the atomic transition frequency lu fixed while tuning the 
cavity resonance u> c . On the other hand, to measure the 
forward emission spectrum, the cavity resonance lo c is 
fixed and the atomic transition frequency ujq is varied. 

The side emission spectrum here is the same as the 
spontaneous emission spectrum calculated elsewhere [HJ . 
The forward emission spectrum is what we expect to 
measure by an ideal detection system at the output of 
the cavity in the forward direction, and has not been 
presented previously, to our knowledge. For zero atom- 
cavity detuning A = 0, where the atom and cavity 
resonances are degenerate, both the side emission and 
forward emission spectra show the normal-mode split- 
tings, which however, are different. The splittings are 
Alo s = 2-^/[<7q + 2g^n(n + 7)] 1 / 2 — k 2 for the side emis- 
sion, and Au>f — 2^/g 2 , — (n 2 + 7 2 )/2 for the forward 
emission, as shown by the thicker red curves in Figs.[3Ja) 
and [3jb) respectively. Both are different from the gener- 
alized Rabi splitting 2g. 

Beyond the energy-splitting difference at zero atom- 
cavity detuning, it is also illuminating to investigate the 
dependence of the energy eigenvalue structure on the 
atom-cavity detuning. Shown in Fig. [3] are plots of the 
spectra, in the strong-coupling regime, for seven differ- 
ent values of atom-cavity detuning A. As |A| increases, 
the vacuum Rabi splitting also increases for both the side 
emission and forward emission spectra. At the same time, 
for the side emission spectra, the cavitylike peak features 
stronger emission and the atomlike peak grows smaller. 
While for the forward emission spectra, however, both 
peaks show the same emission intensity. 



IV. INFLUENCE OF PURE DEPHASING ON 
THE NORMAL-MODE OSCILLATIONS AND 
EMISSION SPECTRA 

Generally speaking, pure dephasing means the decay of 
the dipole coherence without change in the populations 
of the system. Any real transition to other states leads to 
population decay. Thus the pure dephasing is caused by 
virtual processes which start from a relevant state and, 
through some excursion in the intermediate states, return 
to the same initial state. These virtual processes give 
rise to the temporal fluctuations of phases of the wave 
functions, which consequently lead to pure dephasing. 



A. Phase-diffusion model of pure dephasing 

The effects of pure dephasing can be calculated numer- 
ically, based on the Green function formalism by consid- 
ering the microscopic details of various virtual processes 
[33l ] . Instead, for simplicity we treat this problem analyt- 
ically in the phase-diffusion model where the incoherence 
due to elastic collisions or elastic phonon scattering is de- 
scribed by a stochastic model of random frequency mod- 
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FIG. 3: (Color online) Normalized (a) side emission 
spectra, and (b) forward emission spectra of SPS for 
seven different values of atom-cavity detuning A /go = 
(2.4, 1.6, 0.8, 0, -0.8, -1.6, -2.4), given (g , k, 7 )/2tt = 
(8.0, 1.6, 0.32) GHz. 
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ulation, as shown in Fig. [4] replacing the atomic transi- 
tion frequency or the phase of the wave function by an 
instantaneous one 

uj -> U (t) = ujq + f(t) or 

ft (16) 
uj t — > / dt'ui (t') — uJot + ip(t), 
Jo 

where f(t) is the instantaneous deviation of the transition 
frequency due to the elastic collisions or scattering pro- 
cess and <p(t) = /q f(t')dt' is the instantaneous stochastic 
phase of the wave function. 



(|20[) can be solved for the quantum expectation value of 
v(t) exactly [36[. For a Wiener-Levy process, the stochas- 
tic average of the equation satisfies the following differ- 
ential equation 



dt 



(v(t)) = [M - lp Mf\ (v(t)) ■ 



(21) 



The solution to Eq. 
transform form 



(*(*)) 
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exp {zt)N-\z) {v(0)), 
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FIG. 4: (Color online) Schematic diagram of pure dephasing 
process in the phase-diffusion model. 



where the matrix N~ l {z) is the inverse to N(z), which 
itself is given by the formula N(z) = zl — (Mo — jpMf ). 

In Eq. (|22|) . the contour of integration C lies parallel 
to the imaginary axis in the complex z plane, to the 
right of all singularities of the integrands. In order to 
find the time behavior of (v(t)), we have to invert the 
matrix N(z), whose determinant plays an essential role 
because the roots of its secular equation are the poles of 
the integration in Eq. 



We assume the phase of the wave function is a Wiener- 
Levy process [H)]. In this phase-diffusion model f(t) = 
ip(t) is a random, stationary, Gaussian variable with the 
mean value and the meansquare correlation given by 

</(*)> = 0, </(*)/(*')> = 2 lp 5(t-t'). (17) 

The angular brackets indicate a statistical average over 
the random variables of the stochastic process. The 
Markovian nature of the process is reflected by the pres- 
ence of the delta function S(t—t'). The Gaussian property 
is introduced such that all higher correlation functions 
can be obtained from the second-order correlation func- 
tion by permutations and multiplications (3l| . 2^ p is the 
pure dephasing rate. 

Taking into account the pure dephasing modeled by 
the stochastic process, the net change for Eq. ^ is that 
the phase term e lAt should be replaced, such that 

E(t) = -^oe i[At+ ^ (t)I C(i) - 7 .E(t) (18) 
C(t) = ~ig e- t[At+v( - t)] E(t) -nC(t). (19) 

We note that the above equations with stochastic random 
variables are examples of a multiplicative stochastic pro- 
cess, studied intensively by Refs. [H, [35|, Hi| . We solve 
these equations exactly, using the method developed by 
Wodkiewicz [36[ for a multiplicative stochastic process 
described by the following general vector equation 



■if(t)Mi\v(t), 



(20) 



where v(t) is an n-dimensional vector, Mo and Mi are ar- 
bitrary nxn matrices, in general complex and time inde- 
pendent, and f(t) is the random variable of the stochastic 
process described by Eq. (fTT)) . The equations of the type 



B. Normal- mode oscillations and quantum 
efficiency in the presence of pure dephasing 

We proceed to solve the stochastic Eqs. (TTSl and (|19p . 
For simplicity, we consider the case when the emitter and 
the cavity are in resonance, A = 0. The detuning can al- 
ways be put back without difficulty. First, by making the 
substitutions E(t) = E(t)e'^ and C{t) = C(t) e - Kt for 
convenience, we obtain the following simpler equations 



E{t) 



C(t) = ~ iga e rt e- l ^E(t). 



(23) 
(24) 



Then by defining variable Y(t) = e- Tt+l ^C{t) so that 
the differential equation for E{t) does not explicitly de- 
pend on the random variable fit), we obtain a matrix 
equation of the type PD|) for a multiplicative stochastic 
process, with 



W = (y® U|,= 



-I -t ).*-(?! 

(25) 

and a statistically independent initial condition 
(v(0)) T = (E(0),Y(0)). 

The inverse matrix to the matrix N(z) is 



N-\z) 
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z + r + 7 P 
-igo 



z 



(26) 



Plugging Eq. (|26|) back into Eq. ([22|) . using the Laplace 
transform technique and choosing properly the contour of 
integration C, we obtain the quantum expectation value 
of the probability amplitude (E(t)), 
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where g x = ~ ( r + 7 P ) 2 / 4 - 

Similarly, if we define X(t) = e^-^W E(t), while keep- 



ing C(t) unchanged, we obtain 
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where ,g 2 = Vffg = ( r -7p) 2 / 4 - Taking into account the X(0) = £7(0) = E(0) and F(0) = (7(0) = C(0), we 
definitions of E(t) and (7(i), as well as the fact that transform back to E(t) and C(t), 
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The generalized Rabi frequencies for (E(t)) and (C(t)) 
now are different from each other, as compared to Eq. |6]), 
where they were the same for both (E(t)) and (C(t)}. 
This implies the destroying of coherence between the two 
eigenstates of the system, due to the pure dephasing pro- 
cess. We will show this phase-destroying effect on the 
normal-mode oscillations explicitly later in this section. 
We are more interested in finding the influence of the 



pure dephasing on the probabilities \C(t)\ 2 and \E(t)\ 2 , 
or I(t) = \C{t)\ 2 and J(t) = \E(t)\ 2 , because they give 
the normal-mode oscillations and are what one measures 
in experiment. In order to find the equations of motion 
for them, we have to introduce two other one-time func- 
tions H{t) = E{t)C*{t) and H*(t) = E*(t)C(t). The 
equations of motion for these functions are 



J 



H(t) 
d_ \ H*(t) 
dt | I(t) 
J(t) 






V igae-™-^ 



-igoe' 






Tt+iip{t) 



-igoe 
igoe' 



-rt+i(p(t) 
-rt-i<p(t) 



-igoe 1 



ig e-"^-' \ 

a Tt-itp(t) 




(30) 



r 



We still solve these equations assuming the quantum 
emitter is prepared in an excited state E(0) = 1, (7(0) = 
at time to = 0. We solve these one-time functions one 
by one as we did above for solving (E(t)) and (C(t)). For 



example, to find the solution to (I(t)), defining Ui(t) = 
e rt-Mt)H(t), U; = e rt+i ^H*(t), Zi{t) = e 2Tt J{t) and 
keeping I(t) unchanged, we obtain a matrix equation as 
the standard vector form of Eq. (|20p. with 
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( Uj(t) \ 

i(t) 
V zi(t) J 







V 





r 

-Wo igo 

igo -igo 



-ig ig 

igo ~igo 


2r 





1 















1 







I 












V o 








0/ 



(31) 



r 



and a statistically independent initial condition 
(vi(0)) T — (0, 0, 0, 1), where we have used the 
initial conditions at t — 0, as well as the definitions of 
E(t), C(t) and 



In the cavity-QED strong-coupling regime, (4gr§ — 
T 2 ) > r 2 , 7 2 , the solution for (I(t)) is found to be well 
approximated by (see Appendix IB ip 



J 



</(*)> 



- £o_ [T- lp (l+e)/2]t 



7 P (l+3e)*/2 _ Jp 
4ff 



sm(2gi) — cos(2gt) 



(32) 



r 



where e = (r/2g) 2 , and g = ^/g 2 — (r/2) 2 is the general- 
ized Rabi frequency, as defined before. Treating j p /g as 
a perturbation parameter, we kept the order to 0("f p /g) 
in the coefficients and the order to 0(j p e/g) in the ex- 



ponential arguments. 

Similarly, after some tedious algebra, we find the time 
evolutions of (J(i)) and (H(t)) are (see Appendixes IB 21 
and[B3|, 



2# 2 



7p _ g(r - i P /2) 

4ff 



(H(t)) 



!50 p -7p(3+e)t/2 

25 



3r 



7*(l-7e)t/2 



r - 7 P - 



5o 2 

7 P (l-9e)t/2 



sin(2gi) 

r + 2 7p 

25 



1 - 



5o 2 



cos(2gt) 



cos(2gt) + sin(2gi) 



(33) 
(34) 



Finally, the quantum expectation value of the complex 
conjugate of H(t) is just the complex conjugate of its 
quantum expectation value (H*(t)) = (H(t))*. Using 
the definitions of E(t) and C(t), we can easily find the 



,-27* 



(J(t)> and <|C(t)| 2 > 



solutions for (\E(t)\ 2 ) 
e- 2Kt {I{t)). 

Therefore, the probability of finding the system in the 
excited atomic state, including the pure dephasing, is 



(P e {t)) = (\E(t)\ 2 ) 



gQ p -[K+ 7i ,(l+ E )/2]i 

2ff 2 



j7p (l+3e)t/2 



1p s(r - 7 P /2) 



gl 



sin(2gt) 



2.9 



, jC os(2 5 t)^.(35) 
5o 



And the probability of finding the system in the single 
cavity mode with pure dephasing process is 



(P c (t)> = (\C(t)\ 2 ) 



2\ = _£0_ [K+ 7! ,(l+ e )/2]i 



2 5 2 



e 7 P (l+3e) t /2 _ > sin(2 ^) _ cos (2 ff i) 

4 ff 



(36) 



Shown in Figs-JS^a) and[5Jb) are three plots of each prob- 



ability in the presence of pure dephasing, 

The modulation depths of the red-dot and blue- 
triangle curves, with pure dephasing rates 7p/27r = 
(1.0, 2.5) GHz, are reduced, as compared with the black- 
square curves where there is no pure dephasing. The 
normal-mode oscillation frequency seems unaffected be- 
cause we solved for the probabilities only up to first or- 
der in 7 P /g. In fact, it will change slightly from 2g to 
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FIG. 5: (Color online) Probabilities of finding the system (a) in the excited atomic state and (b) in the cavity mode, with 
logarithmic scale for three different pure dephasing rates 7 P /27r = (0, 1.0, 2.5) GHz (black square, red dot and blue triangle 
respectively), given (go, K, 7)/2?r = (8.0, 1.6, 0.32) GHz. 



2<?(1 — 7p/32g 2 ) if we approximate to second order in Consequently, the emission probability of a single pho- 
7 P / g. The normal-mode oscillations are smeared in the ton into the forward beam and the QE with the pure 
presence of the pure dephasing process. dephasing process, as defined before, are 



Kn 2 ( 1 - p-iK-lvt)*- 

(PS)) = 2«/ dt>{\C{t')\*) = ^ \— 



Vqilp) 



K + lv {l + e/2) 



K~ lp e [K + lp (l + e)/2Y + (2gy 



Kg 2 e -[Jr+7p(l+e)/2]t 

!F[K + lp (l + e)/2] 2 + (2gy 



[K + lp (l + e)/2] 



— sm(2gt) + cos(2gt) 
4.9 



2 5 



cos(2gt) - sin(2gt) 



(P oi t -+ oo)) = 4^ (l- (*-7, e )[g + 7p(l + g /2)] 



g 2 K - 7 p e 



[K + lp {l+e)/2Y + (2gf 



(37) 
(38) 



They reduce to our earlier results, given by and below 
Eq. (fTTj) , in the limit j p — > 0. Figure[H{a) are plots of the 
emission probabilities with and without pure dephasing, 
and Fig. [6]Jb) is the QE i] q as a function of the pure 
dephasing rate. The emission probability is also smeared 
for a pure dephasing rate 7 P /27r = 4 GHz compared with 
no pure dephasing. The QE decreases only about 1% as 
the dephasing rate increases from to 4 GHz. 



C. Two-time correlation functions and the 
emission spectra in the presence of pure dephasing 

In order to calculate the emission spectra, we need to 
find the two-time correlation functions because the emis- 
sion spectra in the long-time limit are proportional to 



the Fourier transform of their convolutions. The two- 
time correlation functions are defined as follows: 

Q(t,t') = E(t)E*(t'), R(t,t') = C{t)E*(t') (39) 
F(t,t') ee E(t)C*(t r ), G(t,t') =C(t)C*(t'). (40) 

Of these Q(t,t') and G(t,t') are required for calculating 
the emission spectra. 

The quantum regression theorem [33, [38| , which pro- 
vides a framework to calculate two-time correlation func- 
tions, states that the equations of motion for two-time 
correlation functions Q(t, t') and R(t, t') and F(t, t') and 
G(t, t') with respect to variable t obey the same equations 
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FIG. 6: (Color online) (a) Emission probability without and with pure dephasing rate 7 P /27r = 4 GHZ (black square and red 
dot curves respectively), (b) The QE r\ q as a function of the pure dephasing rate. Other parameters the same as in Fig. [2] 



of motion as those for E(t) and C(t), respectively, 

d t Q(t,t') = -ig e~ rt e^R(t,t'), (41) 

d t R(t,t') = ~tg e rt e'^Q(t,t'), (42) 

d t F(t,t') = -tg e~ rt e^G(t,t'), (43) 

d t G(t,t') = ~tg e rt e~^F(t,t'), (44) 

but now with initial conditions 

Q(t = t\t') - \E(t')\ 2 ^ J(t>), 

R(t = t',t') = C(t')E*(t') = H*(t'), (45) 

F(t = t',t') = E(t')C*{t') = H(t'), 

G(t = t\t') = \C(t')\ 2 = I(t'), (46) 



which are already solved and given explicitly by Eqs. 
to (El. 



We now specialize to the case t > t' and define 
t = t 1 + t. We are particularly interested in the expecta- 
tion values of (Q(t' + r, t')) and {Git! + r, t')) as pointed 
out before. The solutions for their expectation values, 
according to the quantum regression theorem, have the 
same forms as the solutions for the one-time averages of 
(E{t)) and (C(t)) in Eqs. (J27J) and (J2HJ), with the initial 
conditions given above: 



(Q(t' + T,t')) = e -( r +7 P W2 

(G(t' + T,t')) = e {T+ ^ T ' 2 



cos(5ir) + ~ l ~ 7p sin(gir) 
cos(5 2 r) - F lp sin( 52 r) 

2ff2 



(I(t')) 



Wo . , , 
— su%it) 

9i 

sm(.g 2 T) 

92 



(#*(*')} j, (47) 
(H(t')}\. (48) 



By substituting E(t') = E(t')e^' , C(t') = C(t')e Kt ' (C(f + r)C* (t')) . The side emission and forward cmis- 
and the initial conditions into Eqs. (f47|) and (|48p . we sion spectra are then 
obtain the explicit solutions for (E(t ! + t)E* (t')) and 
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27 

— Rc 



dTe iQT 



dt' {E(t' + r)E*(t')) 



7/tt 



[{K + lp )/2-i£l] 2 +g\ 

r - -f P 



-3r 



2 5 2 -[if + 7p (3 + e)/2] (r/2 + 7p ) 



2 (if + 7p + 4 7p e) [if + 7p (3 + £ )/2] 2 + (2. 9 ) 2 if + 2 7p - 4 7p£ j 

2 f r^W,^, _.™ 1 f 1 if- V K / 3o 2 + 7 P [l + (l-2.g 2 / ffo 2 ) e /2] 



55 Re J ( 7 /tt) (re + 7 P ~ ify I I 

9 2 \[(K + lp )/2-*n] 2 +gl \K-J P e 



[if + 7p(l + £)/2] 2 + (2«?) 2 



(49) 



Sps(fi) = — Re<l / dre 



JSlT 



dt' (C(t' + r)C*{t')) 



k/t: 



9 Z ' { [(if + lp )/2~in] 2 +gl 



3r 



2<? 2 -[if + 7p (3 + e)/2] (r/2 + 7p ) 



r - 7 P 



2 (if + 7 P + 4 7p e) [if + lp (3 + e) /2] 2 + (2«?) 2 



if + 2 7p - 47 p e 



ffo Re J (7 + 7p - ifi) 



if + 7 p(l + e/2) 



[(if + lp ) 12 - in] 2 + g\ J 1 K - lp e [ K + lp (1 + e) /2] 2 + {2gf 



(50) 



Shown in Figs.[7Ja) and[7fb) are the side emission and 
forward emission spectra in the long-time limit, for dif- 
ferent dephasing rates 7 p , while other parameters are the 
same as those in Fig. [2] The effect of pure dephasing is 
twofold. The phase fluctuations decrease the peak inten- 
sities of the spectra, and broaden the linewidths of the 
two peaks and hence smear out the splittings, which cor- 
respond to damping rates of the Rabi oscillations in the 
time domain as shown in Fig. [Sj This effect is further 
seen to increase with increasing values of pure dephasing 
rate 7p . 

V. CONCLUSION 

We derived analytical formulas for the side and forward 
(useful cavity output) emission spectra of single-photon 
sources in the cavity-QED strong-coupling regime. We 
also studied the influence of the pure dephasing process 
on the emission spectra and the QE, in the case that 
the pure dephasing rate is significantly less than the co- 
herent coupling rate, that is, up to first order in j p /g. 
These results should be useful in analyzing photolumi- 
nescence spectra from strongly coupled semiconductor- 
QD microcavities, where pure dephasing cannot always 
be assumed negligible because often temperature tuning 
of the QD has to be used to tune through cavity reso- 
nance [24], [2f| . One can use this method, for exam- 
ple, to model the time jitter of solid-state SPS, where 
the excited state of the QD or color center in diamond 
is often populated by spontaneous phonon emission, by 
averaging over nonradiative relaxation time. One may 
also calculate the two-photon interference visibility as- 
suming having two independent but identical SPS and 
investigate how the pure dephasing processes affect the 
indistinguishability of the emitted single photons. 
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APPENDIX A: FORWARD EMISSION AND 
SIDE EMISSION SPECTRA 

For a stationary and ergodic process, the Wiener- 
Khintchine theorem [3l| states that the spectrum is given 
by the Fourier transform of the two-time correlation func- 
tion of the radiated field. One can easily generalize the 
definition of the Wiener-Khintchine spectrum to that of a 
nonstationary spectrum appropriate in this case. We de- 
fine both the side emission, the spontaneous emission of 
the excited emitter into the free-space other than the cav- 
ity (side modes or leak modes), and the forward emission, 
the emission of single photons through the cavity mirror 
into a single wave-packet, outward-traveling wave. We 
recognize from Eq. ((4]) in the text that OM) is propor- 
tional to the Fourier transform of the probability ampli- 
tude C(t'), 

Og(i) - -iBt f dt'e^ k ~^)t' C (t'). (Al) 
Jo 

We define the spectrum as the absolute value squared 
of the Fourier transform of the probability amplitude in 
the long-time limit, which is proportional to the Fourier 
transform of the convolution of the probability ampli- 
tude, as will be shown later. Therefore, the forward emis- 
sion spectrum is given by 

S fe (oj - cu c ) = lim D(w c ) \0^(t)\ 2 , (A2) 
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spectra (a) and forward emission spectra (b) for three val- 
(black square, red dot and blue triangle respectively), given 



where we have changed the probability amplitude from 
Oj:(t) to Ou(t) by using the density of states for the one- 
dimensional photon reservoir D(uj c ) — L/2nc (38|. Using 
the solution to the probability amplitude C(t) and the 
expression of Oz(t) in Eq. (|Aip . we can calculate the 
spectrum 



ssE={2lf^ dt\E(t)\ 2 ^j S SE - (A7) 

where f2' = u> — luq is the side emission frequency centered 
at the atomic transition frequency ujq. 



S fe {uj-uJc) = lim D(oj c )\B(oj c )\ 2 / dt'e^-^' 
xC(t') f dt"e-' l ^- UJ ^ t "c*(t").(A3) 



Then using the definition of the decay rate of the intra- 
cavity field k = ttD(ui c )\B(ui c )\ 2 38] and defining a new 
variable il = lu — lu c , the forward emission frequency cen- 
tered at the cavity resonance uj c , and r = t' — t", we 
obtain the forward emission spectrum 



x Rc 



2 k 



dre 1 



dt'C(t' +T)C*(t') 







The normalized forward emission spectrum is 

-l 



sfe 



2k I dt\C(t)\' 
o 



S. 



FE- 



(A4) 



(A5) 



Similarly the side emission spectrum and the normal- 
ized side emission spectrum, in the long time limit, are 
given by 



SseW) = ^ 

7T 



x Rc 



r°° r r°° 1 ^ 

J dre in ' T J dtE{t + t)E* (t) j 



(A6) 



APPENDIX B: THE APPROXIMATE 
SOLUTIONS FOR THE EXPECTATION VALUES 

OF I(t), J(t) AND H{t) 

1. Approximate roots of the secular equation of 
the matrix Nj(z) and the solution for in the 

limit (4<? 2 -r 2 ) »r 2 , 7 2 

From Eq. (|3Tj) in the text, define matrix M = Mo — 
jpMi, given explicitly by 



r-7 P -ig igo 
-i go i go 

igo -igo o 2r 



(Bl) 



Then the matrix Ni(z) = zl — M and its determinant 
are 



z + 7 P - r 









Wo -Wo 

-igo igo 



igo 




-igo 


igo 


z 








z-2T 



(B2) 



and 



det [Ni{z]] = {z + lp -T) 

[z{z-2Y) (z + 7p -r)+4.g 2 ^-r)] 



(B3) 



The secular equation is given by the vanishing of the 
determinant Eq. (|B3[) . which reduces to a cubic equation, 
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for z\ 



1p, 



(z - T) [(z - Yf + lp (z - r) + Agl - r 2 ] = 7p r 2 ,(B4) 

which is the standard Torrey equation p^. |40| . This cu- 
bic equation can be solved exactly [4l|, although only 
in an implicit form. As Torrey has pointed out, in the 
special case of interest, this equation has a relatively sim- 
ple explicit solution. We solve it in the strong-coupling 
regime, (4<?q — T 2 ) ^> T 2 , 7 2 , in which case there are two 
kinds of roots. The first of these follows the assumption 
that (z — r) 2 is small compared with (4<7q — T 2 ), allow- 
ing one to rearrange the cubic equation (|B4|) and solve 
by iteration 



7 P r 2 



4.9 2 - I* 



(z-T)(z-T + 7p ) 



4 5o 2 - r 2 



Z2 



O 



(B5) 



where e = (T/2g) 2 , and note that g = yj 'g 2 - (T/2) 2 . 
The second kind of root occurs when (z — T) 2 is as large 
as (4<7q — T 2 ), but with opposite sign. The cubic equa- 
tion (IB4I) can be written as 



(z-r) 2 + 4.g 2 -r 2 = - 7p (z-r) 



{z-Tf 



.(B6) 



To first order in 7p , the factor (z — T) 2 on the right-hand 
side, Eq. (|B6[) can be replaced by -(4g 2 -T 2 ). This gives 



a quadratic equation for (z — T), (z — T) 2 + 7p (l + e)(z — 
r) + 4g 2 = 0, whose solutions are the third and fourth 
roots 



z 3A » r - ^(1 + e) ± i2g^ 1 - ( 7p /4g) 2 + O 



T-^(l + e)±i2g + 



(B7) 



The inverse matrix to the matrix Nj(z) in Eq. (|B2[) is 



^(z) 



det [A7(z)] I "31 "32 "33 ™34 



with 



"3i = -igo ( z - 2r) (z + 7 P - r) 

n 3 2 = «ffo - 2r) (z + 7 P - r) 

«33 = {z + lp -T)[{z-2T){z + lp -T) + 2g 2 } 

"34 = 2g 2 (z + 7 P - r) 

where we only calculate the elements of the third row 
of Nf (z) because they are required to calculate (I(t)), 
which is then 



(I(t)) = 



c 

A 

2g 2 



dz 
~2~ki 



"34 



(z - zi) (z - z 2 ) (z - z 3 ) (z - z±) 



dz 
2ni 



25o 2 



(z - z 2 ) (z - z 3 ) (z - Zi) 



,[r-7p(l+e)/2]* 



45 



sm{2g£) — cos(2gt) 



(B8) 



where we have used the initial condition that (vi(0)) T = 
(0, 0, 0, 1) Treating j p /g as a perturbation parameter, 
we kept the order to 0(j p /g) in the coefficients and the 
order to 0{^f p e/g) in the exponential arguments. 



2. Approximate roots of the secular equation of 
the matrix Nj(z) and the solution for (J(t)) in the 

limit (4 5o 2 -r 2 ) »r 2 , 7 2 



As it is clear from the definition of vi(t), we can only 
obtain the solution for (I(t)) in the above calculation. In 
order to obtain the solution for (J(t)), we have to derive 
another equation of the type Eq. (|2T)|) with the following 



definitions of the vector and matrices: 



Mo = 



/ Uj(t) 
Wj(t) 

V m 















|.M,= 




i 



























/ 




(B9) 



and with the initial condition (v,j(0)) T — (0, 0, 0, 1), 
where Uj(t) = e - r *-^(*) H (t) , U}(t) = e- rt+i ^H*(t), 



Wj{t) 



-2Tt 



I(t). 



The calculation of Nj(z) is almost the same as the 
calculation in Appendix IB II The matrix M = Mn — 
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IpMl is 



The inverse matrix to the matrix Nj(z) is therefore 



M 



-r - 7 P -ig «.go 

-z.go igo -2r 
iffo ~«5o 



(BIO) 



Then the matrix Nj(z) = zi — M and its determinant 
are, respectively, 



det [iVj(z)] 



• • • 

• • • 

• • • 



7141 H42 "43 U44 



/ z + 7 P + r ig -ig 

2 + 7 P + T -ic/o «5o 

igo -igo z + 2T 

\ -ig ig z 



with 



and 



det [Nj(z)} = (z + 7p + r) 

[z (z + 2r) (z + 7p + T) + 4.g2 (z + T)] 



(Bll) 



which is the same as Eq. (|B3I) provided that we change T 
to — r. So the roots of the secular equation of the matrix 
Nj{z) are 



"41 = igo (z + 2T) (z + 7 P + r) , 

"42 = -igo (z + 2T) (z + 7 P + r) , 

"43 = 2gQ (z + 7 p + T) , 

"44 = (z + 7 P + T) [(z + 2T) (z + 7 p + T) + 2gl] 



Zx » -r - 7p, z 2 -r + 7 P £ + o 



z 3 , 4 ^-r-^(l + £ )±z2,g + 



(B12) 



where we only calculate the elements of the fourth row 
of Nj(z) because they are required to calculate (J(t)), 
which is then 



dz 



"44 



c 2-Ki (z - zi) (z - z 2 ) (z - z 3 ) (z - z 4 ) 



dz cZt (z + 2r) (z + r + 7p ) + 2gg 

C 27TJ (z - z 2 ) (z - z 3 ) (z - Zi) 



gQ e -[r+ 7p (l+ E )/2]t ^ e 7 P (l+3e)t/2 

2# 2 



> _ 5 (r - 7p/2) 
4ff 



3o 2 



sin(2<?£) 



1 



^- ) ccsl 2.7/1 



(B13) 



where we have used the initial condition that (vj(0)) = 
(0, 0, 0, 1) and kept the order to O {j p /g) and O (T/g) in 
the coefficients and the order to O (j p e/g) in the expo- 
nential arguments. 



definitions of the vector and matrices: 



v H {t) 



Mo = 



( H(t) 
U H {t) 
W H (t) 

V z H {t) 





/ 
2 
10 

\o o o i 



3. Approximate roots of the secular equation of 
the matrix Nh(z) and the solution for (H(t)) in the 



limit (Agl - r 2 ) > r 2 



7 P 



In order to obtain the solution for (H(t)}, we derive 
another equation of the type Eq. (f2"0")) with the following 



and with the initial condition (v H (0)) J = (0,0,0,1), 

where U H (t) = e i2 ^H*(t), W H {t) = e - rt+lv ^I{t), 
Z H (t) = e rt+i ^J{t), 

The matrix M for the vector vu(t) is 



M = M Q - lp Ml = 









-igo 


igo 





-47 P 


igo 


-igo 


-igo 


igo 


-r-7p 





igo 


-igo 





r -7p 
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Therefore 



N H {z) 



( z ig -ig 

z + 4 7p ~ig ig 

igo -igo z + (r + 7p ) o 

V ~igo wo o z - (r - 7p ) 



and the determinant of the matrix iVy (z) is 

det[A^(z)] = z(z + 4 7p ) [(z + 7p ) 2 - T 2 

+ 4. 9o 2 (z + 7p )(z + 2 7p ). (B14) 

The secular equation is given by setting det [JVh (z)] = 0, 
which gives 



The second kind of root occurs if (z + 7p ) (z + 2 7p ) is as 
large as {Ag^ — T 2 ), but has the opposite sign. Then the 
alternative rearrangement of Eq. (|B15j) is 



(z + 7p )(z + 2 7p ) + (4 5o 2 -r 2 - 4 7 2 ) 

z - 2 7p 



(B19) 



(z + 7p ) (z + 2 7p )_ 

To first order in 7p the factor (z + 7p ) (z + 2 7p ) on 
the right hand side of Eq. (|B20p can be replaced by 
— (4g 2 — T 2 ). This gives a simple quadratic equation for 
z, z 2 + 7p (3 + e)z + 4g 2 - 4 7 2 + 2 7 2 (1 - e) = 0, whose 
solutions arc the third and fourth roots 



(z + 7p ) (z + 2 7p ) + 4 ffo 2 (z + 7p ) (z + 2 7p ) - 

T 2 (z + 2 lp f - 4 7p 2 (z + lp f = -4 7p T 2 . (B15) 

In the most general case, no simple factorizations oc- 
cur, and a quartic equation must be solved. Again the 
roots are implicit in the general case [42j], but explicit 
in the strong-coupling regime. Similarly, there are two 
kinds of roots in the strong-coupling regime. The first 
of these follows from the assumption that both (z + 7p ) 2 
and (z + 2 7p ) 2 are small compared with (4gp — T 2 ) , in 
which case it is natural to rearrange Eq. (|B15|) into the 
form 



(z + 7p ) (z + 2 7p ) 
x [(z + 7p ) (z + 2 7p ) + (4 3o 2 - T 2 - 4 7p 2 )] 



-4 7p 2 r 2 , 

(B16) 



where we used the assumptions (4<7q — T 2 ) » T 2 , 7p , 



z+7p 



1 and 



z+7p 



(z + 7p ) (z + 2 7p ) 

-4 7p 2 r 2 

4 ff 2 - P - 4 7p 2 



1 



1. Then 



(z + 7p) (z + 2 7p ) 



r 2 -4 7 2 



(B17) 



which is solved by iteration. The roots are 

3" 



zi « - 7p (l + 4e) + 



z 2 » -2 7p (l - 2e) + O 



(3 + s)± l 2g^l-( lp /g) 2 + 



-3,4 * -y 

« -y (3 + e)±i2g + 



The inverse matrix to Nh(z) is 



det[iV(z)] 



/ nn "12 "13 "14 

• • • • 

• • • • 

• • • • 



with 



(B20) 



nn 
"12 

"13 
"14 



(z + 4 7p ) (z + 7p ) 2 -r 2 



+ 2ff 2 (z + j p ) 



= o, 



-ig (z + 4 7p ) (z + 7p - T) . 

igo {z + 4 7p ) (z + 7p + r) . 



(B18) Therefore is given by 



(H(t)) 



Jc 



dz 
~2m 



,zt 



ifO 7p(3 +e)t/2 
2/ 



"14 

(z - zi) (z - z 2 ) (z - z 3 ) (z - Zi) 

r ^£ e 7p(l-7e)t/2 _ L_ 

tg g 



igo (z + 4 7p ) (z + 7p + r) 

) c (z - Zi) (z - z 2 ) (z - z 3 ) (z - z A ) 

> e -7p(l~9e)t/2 r + 2 7p 



dz 
2m 



2g 



- cos(2gt) 



- z s) [z 
- sm(2gt) 



(B21) 
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where we have used the initial condition that (vh(0)) = in the coefficients and the order to O {jps/g) in the ex- 
(0, 0, 0, 1) and kept the order to O (j p /g) and O (T/g) ponential arguments. 
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